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Abstract 

Metal nanowires exhibit a number of interesting properties: their electrical conduc- 
tance is quantized, their shot-noise is suppressed by the Pauli principle, and they are 
remarkably strong and stable. We show that many of these properties can be under- 
stood quantitatively using a nanoscale generalization of the free-electron model. Possible 
technological applications of nanowires are also discussed. 

Introduction Metal nanowires represent nature's ultimate limit of conductors down to a 
single atom in thickness. In the past eight years, experimental research on metal nanowires 
has burgeoned [1-13]. The simplest model of a metal is the free-electron model H], which 
already describes many bulk properties of simple monovalent metals semiquantitatively. In 
this article, we discuss our generalization of the free-electron model to describe nanoscale 
conductors [15-22]. 

A remarkable feature of metal nanowires is the fact that they are stable at all. Fig. |l| 
shows electron micrographs by Kondo and Takayanagi JsJ illustrating the formation of a gold 
nanowire. Under electron beam irradiation, the wire becomes ever thinner, until it is but 
four atoms in diameter. Almost all of the atoms are at the surface, with small coordination 
numbers. The surface energy of such a structure is enormous, yet it is observed to form 
spontaneously, and to persist almost indefinitely. Even wires one atom thick are found to 
be remarkably stable |j, ^j, [l^] . Naively, such structures might be expected to break apart 
into clusters due to surface tension [^3[, but we find that electron-shell effects can stabilize 
arbitrarily long nanowires p2| . 

A crucial clue to understanding the physics of metal nanowires is the observed correlation 
between their electrical and mechanical properties. In a seminal experiment ^ carried out 
in 1995, Rubio, Agra'ft and Vieira simultaneously measured the electrical conductance and 
cohesive force of an atomic-scale gold wire as it formed and ruptured (see Fig. ^, left panel) . 
They observed steps of order Go = 2e 2 /h in the conductance, which were synchronized with a 
sawtooth structure with an amplitude of order InN in the force. Similar results were obtained 
independently by Stalder and Diirig Q . Note that the tensile strength of the nanowire in the 
final stages before rupture exceeds that of macroscopic gold by a factor of 20, and is of the 
same order of magnitude as the theoretical value in the absence of dislocations ||. This is 
consistent with the recent finding of Rodrigues, Fuhrer, and Ugarte that such nanowires are, 
in fact, typically free of defects in their central region J13| . 

The standard description of nanoscale cohesion, pioneered by Landman and coworkers , 
is via molecular dynamics simulations which utilize short-ranged interatomic 

potentials suitable to describe the bulk properties of metals. However, such an approach 
appears problematic when applied to metal nanowires, in which electron-shell effects |nj 
due to the transverse confinement are likely to be important. On the other hand, atomistic 
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Figure 1: Transmission electron micrographs showing the formation of a gold nanowirc |J 
(image courtesy of Y. Kondo): (a) an image of Au(OOl) film with closely spaced nanoholes, 
the initial stage of the nanowire; (b) a nanowire four atoms in diameter, resulting from 
further electron-beam irradiation. 

quantum calculations using, e.g., the local-density approximation, are restricted to such 
small systems that their results can not really be disentangled from finite-size effects pQ ]. 
An alternative approach, developed by our group, is to replace the discrete ionic coordinates 
by a coarse-grained jellium background, in order to be able to treat the electronic degrees of 
freedom correctly. We have argued ]I3 that an atomic-scale contact between two pieces of 
metal can be thought of as a waveguide for conduction electrons (which are responsible for 
both electrical conduction and cohesion in simple metals) : Each quantized mode transmitted 
through the contact contributes 2e 2 /h to its conductance and a force of order sf/Xf (roughly 
InN) to its cohesion, where Xf is the de Broglie wavelength of an electron at the Fermi energy 
£f (see Fig. ^, right panel). To my knowledge, our approach is the only one in which the 
observed correlations between the cohesive and conducting properties of metal nanowires 
have been explained within a single theoretical model. 

The paper is organized as follows: The free-electron model of nanoscale conductors is 
introduced in the next section, followed by a discussion of quantum transport, including the 
effect of realistic contacts to the nanowire. Nanoscale cohesion is then analyzed within our 
model, followed by a discussion of the remarkable stability of nanowires. The paper concludes 
with some comments about the technological promise of metal nanowires. 

Free-electron model We investigate the simplest possible model Jl^, [l6| for a metal 
nanowire: a free (conduction) electron gas confined within the wire by Dirichlet boundary 
conditions. A nanowire is an open quantum system, and so is treated most naturally in terms 
of the electronic scattering matrix S. The Landauer formula 0] expressing the electrical 
conductance in terms of the submatrix S12 describing transmission through the wire is 

G =¥/ deZ l^ Tr {^ (e)5i2(e) } T = ¥5: r " ( ^ ) ' (i) 

where f(e) is the Fermi-Dirac distribution function and the transmission probabilities {T n } 
are the eigenvalues of S\ 2 Si2- The conductance of a metal nanocontact was calculated exactly 
in this model by Torres et al. j29| . The appropriate thermodynamic potential to describe the 
energetics of such an open system is the grand canonical potential Q: 

fi = -~ J deg{e) In (l + c -^-">) T =° J^deg(e)(e - e F ), (2) 
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Figure 2: Left: Measured electrical conductance and cohesive force of a gold nanowire 
H] (image courtesy of N. Agrai't). Right: Calculated conductance and force of a metal 
nanowire, modeled as a constriction in a free-electron gas with hard walls ]l5| . Note that 
ep/Xp ~ 1.7nN in gold. 



where /3 is the inverse temperature, [i is the chemical potential of electrons injected into the 
nanowire from the macroscopic electrodes, and g(e) is the electronic density of states (DOS) 
of the nanowire. The DOS of an open system may be expressed in terms of the scattering 
matrix as J3(| 

^) = ^{^( £ )f -H.c). (3) 

This formula is also known as the Wigner delay. Thus, once the electronic scattering problem 
for the nanowire is solved, both transport and energetic quantities can be readily calculated 
[ |l5| , |l7) . Electron-electron interactions can be included at the mean-field level in this 
model in a straightforward way Jl(], |l^, E3|, but do not alter our main conclusions. 

Quantum Transport Evaluating the transmission probabilities {T n } in the WKB approx- 
imation for an axially-symmetric nanowire |15[ , the conductance calculated from Eq. (|j]) is 
shown in the upper-right panel of Fig. ^| Plateaus in the conductance at integer multiples of 
Go are evident, with some rounding of the steps due to tunneling. Some integers are absent, 
reflecting the degeneracies associated with axial symmetry ^] . 

Conductance steps of size Go were first observed in quantum point contacts (QPCs) fabri- 
cated in semiconductor heterostructures and are a rather universal phenomenon in metal 
nanowires [1-4], even being found in contacts formed in liquid metals Q. The precision of 
conductance quantization in metal nanowires is poorer than that in semiconductor QPCs due 
to their inherently rough structure on the scale of the Fermi wavelength Xf, which causes 
backscattering [|17| , and due to the imperfect hybridization of the atomic orbitals in the con- 
tact, especially for multivalent atoms M. For this reason, a statistical analysis of data for a 
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Figure 3: (a) £alculated conductance histogram [171 ; (b) calc ula ted mean shot noise (sj) 
(grey squares mM), together with experimental datafrom Ref. |l0[ (black circles); (c) mean 
transmission probabilities (T n ) |18| . The error bars indicate the standard deviations of the 
numerical results over the ensemble and the experimental errors, respectively. 



large number of contacts is often made ^, ^, ^, [uj, jnj, resulting in a conductance histogram 
[see Fig. |(a)]. 

To model quantum transport in gold nanowires, where there are no "missing integers" 
in the conductance histogram [J 0, jlQ], geometries without axial symmetry were chosen, 
and weak disorder, corresponding to a mean-free path }zfI — 270, was included both in 
the nanowire and in the electrodes neighboring it [p_7| . The transmission probabilities were 
calculated by solving Schrodinger's equation using a recursive Green's function algorithm 
P7[ . Averaging over different contact shapes and impurity configurations, we obtained the 
histogram shown in Fig. |^(a), which is very similar to typical experimental histograms for 
gold Jl], [], |l(|. The effect of disorder is twofold the conductance peaks are shifted 

downward due to backscattering, and the peaks are broadened due to universal conductance 
fluctuations, filtered by the nanowire. 

Recently, additional information on quantum transport in metal nanowires has been ob- 
tained from experiments on shot noise [jlo|. Shot noise is the term used to describe the 
temporal fluctuations of electric current arising from the discreteness of the electric charge 
e. In 1918, Schottky showed that if the arrival times of charge carriers are uncorrelated, the 
shot-noise spectral power Pi — 2el, where / is the time-average current. However, in a quan- 
tum conductor with a finite number of transmitted modes, the shot noise is suppressed below 
the Schottky value due to anticorrelations induced by Fermi-Dirac statistics. The suppression 
factor at zero temperature is given by |l^] 

Sl = fL = ^Ml-T n ) 

Sl ~2el- £„T n • W 

Fig. ^b) shows the measured shot noise (solid circles |HJ) of gold nanowires as a function of 
their conductance. The pronounced suppression of si for wires with conductances near integer 
multiples of Go reveals unambiguously the quantized nature of the electronic transport. We 
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computed |1S| the mean and standard deviation of si and T n as functions of G (grey squares 
in Fig. [}]) from the numerical data used to generate the conductance histogram in Fig. [](a). 
The agreement of the experimental results for particular contacts and the calculated distri- 
bution of si shown in Fig. ^(b) is extremely good: 67% of the experimental points lie within 
one standard deviation of (si) and 89% lie within two standard deviations. It should be em- 
phasized that no attempt has been made to fit the shot-noise data; the numerical data of Ref. 
Jl7[ , where the length of the contact and the strength of the disorder were chosen to model 
experimental conductance histograms for gold, have simply been reanalyzed to calculate (sj). 
The 97% suppression of shot noise for nanowires with a single quantum of conductance (i.e., 
wires one atom thick) suggests that such wires could be useful for low-temperature/low-noise 
applications, such as quantum computing. 

Metallic nanocohesion The cohesive force of the nanowire is F — —d£l/dL, where L 
is the length of the nanowire. We assume that the volume per atom is conserved under 
elongation (ideal plastic deformation), so that the deformation occurs at constant volume 
(for alternative constraints, see Refs. jl9"| , ^l|). While the conductance is determined by 
the transmission probabilities, Eqs. (^) and(g|) indicate that the energetics of a nanowire 
are determined by the scattering phase shifts. Evaluating the phase shifts in the WKB 
approximation, performing the energy integral in Eq. (^|) at T = 0, and taking the derivative 
with respect to elongation Jl^] , one finds the force shown in the lower-right panel of Fig. ^. 
The correlations between the force and conductance are striking: as the wire is elongated 
and its diameter decreases, increases along a conductance plateau, but decreases sharply 
when the conductance drops. Each transmitted mode acts like a delocalized metallic bond, 
which can be stretched and broken. 

The calculated force is remarkably similar, both quantitatively and qualitatively, to the 
measured force for gold nanowires, shown in the lower-left panel of Fig. [| Inserting the 
value bf/Af — 1.7nN for gold, we see that both the overall scale of the force for a given 
value of the conductance and the heights of the last two force oscillations are in quantitative 
agreement with the experimental data. One discrepancy is that the jumps in both force and 
conductance are less abrupt than in the experimental curves, possibly because we considered 
only geometries that change continuously with elongation. 

In order to separate out the mesoscopic sawtooth structure in the force, associated with 
the opening of individual conductance channels, from the overall (macroscopic) trend of the 
contact to become stronger as its diameter increases, it is useful to perform a systematic 
semiclassical expansion |l| H of the DOS, 3(e) = g(s) + 5g{e), where q is a smooth average 
term, referred to as the Weyl contribution, and 8g(e) is an oscillatory term, whose average is 
zero. For the free electron model with Dirichlet boundary conditions, the Weyl term is fe2] 




where k = \^2me/h, V is the volume of the wire, A its surface area, and C the integrated mean 
curvature of its surface. The oscillatory contribution Sg(e) to the DOS may be approximated 
as a Feynman sum over classical periodic orbits a la Gutzwiller fell, p2{: 



where S v is the classical action of a periodic orbit, Q v is a phase shift determined by the 
singular points along the classical trajectory, and A v is an amplitude depending on the 
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Figure 4: Density of states g(ep) of a cylindrical wire (upper diagram) and zero temperature 
stability coefficient a (lower diagram) versus the radius Ro of the unperturbed wire. The 
wavevector of the perturbation is qRo = 1. Dashed curve: Weyl contribution to a. 



stability, symmetry, and period of the orbit. Using g(e) in Eq. (|2|), one can derive a Sharvin- 
like formula for the force 

The first term in F is the surface tension. The oscillatory mesoscopic correction SF may be 
calculated with the aid of Eq. Under reasonable assumptions about the geometry, it can 
be shown J^] that the amplitude of the force oscillations is universal: 

rmsSF = 0.58621 e f /\f. (8) 

In more realistic models including electron-electron interactions jl6|, [l9| |2lj] and self- 
consistent confining potentials pjjfl , the surface tension is typically reduced compared to 
Eq. ([?]), but the force oscillations are essentially the same as in the free-electron model. 

Stability of Nanowires A cylindrical body longer than its circumference is unstable to 
breakup under surface tension [E3|, a phenomenon known as the Rayleigh instability. How 
then to explain the durability of long gold nanowires [c.f. Fig. [j](b)], the thinnest of which 
have been shown to be almost perfectly cylindrical in shape? The key is the quantum 
corrections [^2[ to the classical stability coefficients. 

Only axially-symmetric deformations can lower the surface energy of a cylindrical object, 
and thus lead to an instability (23) . Any such deformation may be written as a Fourier series 

R(z) = R + [ X dqb{q)e lqz , (9) 
J —00 

where R(z) is the radius of the wire at z, Ro is the radius of the unperturbed cylinder, and 
b(q) is a complex perturbation coefficient. Using Eqs. (^) and ^ in Eq. (^|), one obtains the 
following expansion [ p2[ 

r 00 

Q[b] = Q[0] + / dqa{q)\b{q)\ 2 + 0(b 3 ), (10) 
Jo 
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Figure 5: Stability of metal nanowires as a function of radius and temperature. Shaded 
regions indicate stability with respect to small perturbations; unshaded regions denote 
unstable configurations. Here Tp is the Fermi temperature and Ro is the mean radius of 
the wire. The quantized conductance values of the stable wires are indicated. 



where the stability coefficient a(q) depends implicitly on Ro and T. The change in the grand 
canonical potential is of second order in b and contributions from deformations with different 
q decouple. If a(q) is negative for any value of q, then fi decreases under the deformation 
and the wire is unstable. 

Let us first discuss the stability of a nanowire at zero temperature. Fig. ^| shows the stabil- 
ity coefficient (lower diagram) and DOS (upper diagram) at the classical stability threshold 
qRo = 1 as a function of Ro. For a straight wire, the transverse motion is quantized, and 
the DOS consists of a sequence of sharp peaks associated with the opening of each succes- 
sive subband. a has sharp negative peaks — indicating strong instabilities — at the subband 
thresholds, where the density of states is sharply peaked. Under surface tension and curva- 
ture energy alone (dashed curve in Fig. Q), the wire would be slightly unstable at the critical 
wavevector qRo = 1, since the curvature term is negative. However, the quantum correction 
is positive in the regions between the thresholds to open new subbands, thus stabilizing the 
wire. Since the oscillatory contribution to a is independent of q, we find that regions of sta- 
bility persist for arbitrarilly long wavelength perturbations, indicating that an infinitely long 
cylindrical wire is a true metastable state if the radius lies in one of the windows of stability. 

With these results, we can construct a stability diagram for metal nanowires (see Fig. []) . In 
the semiclassical approximation, dct/dq > always, so the stability of the wire is determined 
by the sign of a(q — 0). a(0) is a function of two dimensionless parameters, kpRo and 
T/Tf- Regions where «(0) > are shaded dark in Fig. ^, while regions where q(0) < are 
unshaded. The stable regions persist up to extremely high temperatures for several quantized 
conductance values (recall that Tf > 10 4 K for metals), indicating that electron- shell effects 
may stabilize nanowires even for temperatures well above the bulk melting temperature. It 
is important to point out that if a more realistic value of the surface tension jl^, Q 
were used, the stability boundaries would be pushed to even higher temperatures. Thus the 
astounding stability properties shown in Fig. ^ are a very robust prediction of the jellium 
model. 



Conclusions Wires formed from chains of individual metal atoms have a number of prop- 
erties which make them promising for nanotechnology: They are very strong, able to support 
tensions up to ef/^f ~ InN. Contrary to naive expectations, they are extremely stable, de- 
spite their large surface to volume ratio. They are nearly-ideal one-dimensional conductors, 
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and exhibit dramatically-reduced shot noise. One potential application of metal nanowires 
is for integrated-circuit interconnects, due to their high conductance and structural robust- 
ness. The quantum suppression of shot noise in nanowires may also make them useful for 
low-temperature/low-noise applications, such as quantum computing. 
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